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DYNAMICS OF ARTICULATED STRUCTURES 

For purposes of this presentation, an articulated structure is defined as 
an assembly of flexible bodies that may be coupled by kinematic connections 
and force elements that permit large relative displacement and rotation. 
Kinematics of such systems Is defined using one reference frame for each body 
in the system and deformation modal coordinates that define displacement 
fields within flexible bodies. Deformation kinematics are defined by both 
elastic vibration and static correction deformation modes. Linear eJastic 
deformation is presumed; i.e., a linear stress-strain relation is valid and 
relative displacements within each elastic component are small enough so that 
the theory of linear elasticity applies. Coupling of reference and modal 
coordinates leads to a system of nonlinear equations of motion. Methods of 
automatically generating and solving these equations of motion are outlined. 


• Large Displacement and Rotations of Body Reference Frames 
(reference coordinates) 

• Elastic Vibration and Static Correction Deformation Modes 
(modal coordinates) 

• Coupled Nonlinear Equations in Reference and Modal Coordinates 

• Automated Equation Generation and Solution 
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EULER’S THEOREM 


Euler’s theorem guarantees existence of a unit vector u about which an 
x’-y’-z’ reference frame may be rotated by an angle x to bring it from a 
reference x-y-z frame to a general orientation. Components of the unit vector 
u and angle x °f rotation are used to define orientation of a reference frame 
in space. 


EULER PARAMETERS 


A set of four Euler parameters Is defined as p, as shown on the chart. 
These four parameters are not independent, since the vector p must be a unit 
vector in R . The direction cosine transformation matrix from the x’-y’-z' 
reference frame to the x-y-z frame is defined as shown. The quadratic nature 
of terms in the transformation matrix, as functions of Euler parameters, leads 
to attractive properties when writing velocity and acceleration equations that 
are needed in the equations of motion. Furthermore, use of Euler parameters 
avoids singular orientation difficulties that are associated with a set of 
three rotation parameters, such as Euler angles or Bryant angles [1,2]. 
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LUMPED MASS FLEXIBLE BODY MODEL 


A lumped mass finite-element formulation is used to carry out vibration 
and static correction mode analysis of each deformable body in an articulated 
structure. A typical point P* is defined in the undeformed state of the body 
by a constant vector Tq*. During the process of deformation, this point 
undergoes displacement u* in the body reference frame, as shown. Lumped 
masses m^ at each note in the finite-element model are used in defining 
kinetic properties of the flexible body [3,4]. 
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KINEMATICS 


A vector u of modal displacements is defined as a linear combination of a 
set of deformation modes 4> , j = l,***,m. The displacement vector u 1 of point i 
in the body is provided by a projection matrix P*. As noted earlier, the 
direction cosine transformation matrix for the reference frame associated with 
the body is a function of the Euler parameters of that reference frame • 

Finally, the global position vector R 1 of point i on the body is given as 
shown. 


u » a.^ - nodal displacement relative to reference frame 

<j^ ssjssl 9 # • • jtn - deformation modes 

u* = P*u ~ elastic displacment of mode i 

A(p) - direction cosine matrix of reference frame 

p - Euler parameters of reference frame 

R* * R + A(r* + ) 
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VIBRATION AND STATIC CORRECTION MODES 


Boundary conditions must be selected for characterizing deformation of 
flexible components. Since kinematic constraints on bodies in an articulated 
structure often lead to statically indeterminant sets of boundary conditions, 
a statically determinant or underdetermined set of boundary conditions is 
selected for use in vibration analysis. Unit loads associated with deleted 
kinematic constraints are used to define static correction modes [3]. These 
calculations are carried out with any standard finite-element code that is 
capable of generating lumped mass information. Constants that will appear 
subsequently in the equations of motion are calculated using information 
generated within the finite-element code. 


• Select Boundary Conditions for Flexible Components 

• Calculate Natural Vibration Modes 

• Calculate Static correction Modes for Deleted Constraints 

• Calculate Constants for Equations of Motion 
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KINEMATICS (continued) 


The position relationship derived earlier is differentiated to obtain the 
global velocity vector of node i in the body, as shown. The time derivative 
of the direction cosine orientation matrix yields an expression in the time 
derivative of Euler parameters as shown [2,3]. The velocity vector may thus 
be written in matrix form for use in derivation of the equations of motion of 
the system. 
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KINEMATIC CONSTRAINTS 

A variety of kinematic couplings between flexible bodies is derived in 
Refs. 2 and 3. Joint definition £. Pj Cj reference frames are fixed to the 

deformable body to define information required to write kinematic constraint 
equations associated with each joint in the system. Shown on the chart below 
are spherical, revolute, and universal joints, for which constraint 
equations may be found in Ref. 3. 
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KINETICS 


The kinetic energy of a flexible body may be written in terms of time 
derivatives of reference frame generalized coordinates and modal coordinates, 
as shown. Since vectors r* - and matrices E, G, and A are nonlinear functions 
of Euler parameters, the mass matrix of the flexible body is a nonlinear 
function of generalized coordinates, as a result of geometric nonlinearities 
in the system kinematics. 
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FORCES 


The strain energy U of the flexible body may be written explicitly in 
terms of a modal stiffness matrix K ga , as shown. Generalized forces defined 
through direct application of the definition of virtual work lead to nonlinear 
algebraic expressions in generalized coordinates, as shown. These forces 
include both externally applied forces and forces of interaction due to 
compliant couplings between bodies and feedback control actuators. 
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EQUATIONS OF MOTION OF A SINGLE BODY 


The equations of motion of an unconstrained individual flexible body are 
shown below [3] • The system of equations for an articulated structure that is 
made up of multiple bodies connected by kinematic constraints is developed [3] 
using the Lagrange multiplier form of multi-body system dynamics [1,2]. 
Evaluation of individual terms appearing in the coefficient matrix of 
accelerations and on the right side of the equations of motion is derived by 
expanding the expressions shown and calculating constant coefficients 
associated with deformation modes and mass distribution. 
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CONSTANTS FROM FINITE-ELEMENT MODEL 


Nine sets of constant vectors and matrices shown are calculated, using 
data generated in the finite-element deformation analysis of each flexible 
body. These constants are computed using an intermediate processing program 
[4]. 
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NONLINEAR TERMS IN EQUATIONS OF MOTION 


Three typical nonlinear terras appearing in the equations of motion 
presented earlier are shown here, evaluated as linear and quadratic 
expressions in generalized coordinates. All such terms are coded in a 
f lexible-body module of the Dynamic Analysis and Design System (DADS) computer 
code. These terms are evaluated at every time step in numerical integration 
of the coupled system of nonlinear equations of motion. 
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NUMERICAL SOLUTION 


A variable order, variable step size numerical integration algorithm is 
used to compute the solution of differential-algebraic equations of motion for 
articulated structures. Since step size and order selected by the algorithm 
reflect the error tolerance required and the frequency of oscillation that 
develops, integration cost is influenced by selection of deformation modes in 
the model. Numerical results accumulated to date [4,5] show that prudent 
selection of a combination of vibration and static correction modes gives 
reasonable results. Substantial work remains to be done in rational selection 
of these deformation modes. 


• Variable Order, Variable Step Size Numerical Integration 

• Integration Cost is a Function of Frequency Content 


• Mixed Vibration and Static Correction Modes Give Best Results 


FLEXIBLE HINGE DOOR EXAMPLE 


The flexible door structure shown is kinematically coupled to a body that 
is taken to be rigid ground. The revolute joints shown are misaligned so 
that there is no deformation when the door structure lies in the Y-Z plane. 

Any rotation of the door structure leads to deformation of the beam and plate 
structure of the door, which tends to bring it back to the undeformed state. 
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FINITE-ELEMENT MODEL OF DOOR STRUCTURE 

A modest finite-element model of the door structure, using plate and beam 
elements, is shown. 


N20 NI5 NIO 



NI6 Nil N 6 


N * * : Node Number * * 

Plate ; 12 (Membrane + Bending) Elements 

E = 2.0 * I0 12 N/m 2 

Beam ; 10 Beam Elements at each beam 
E» 2.0xl0 n N/m 2 
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STATIC CORRECTION MODES 


Boundary conditions for finite-element analysis are selected so that the 
center point of the bottom hinge is fixed in space and x- and y-coordinates of 
the top hinge point are likewise fixed. Five kinematic constraints are thus 
suppressed, two rotations at the bottom hinge and two rotations and one 
translation at the top hinge. Unit torques and a unit force are applied to 
calculate five static correction modes to represent deformation of the 
structure. Vibration modes are likewise calculated [4]. 



507 



FIRST APPROXIMATE SOLUTION 


Two models were used in a preliminary analysis. The first consisted of 
only five normal vibration modes, with numerical results for the X-coordinate 
of the center of the door shown as a solid line. A five static correction 
mode approximate solution is shown with a dotted line, reflecting much lower 
frequency of vibration of the door structure. To evaluate reasonables of 
predictions, additional modes and combinations of modes are selected. 


— •• 5N SOLUTION 


« 5S SOLUTION 
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SECOND APPROXIMATE SOLUTION 


The nine normal mode solution shown reflects a somewhat lower overall 
vibration frequency, but is still unreasonable. As shown by the dotted and 
dashed curves, adding four normal vibration modes to the five static 
correction modes yields only a slight change in the prediction obtained from 
purely five static correction modes. This suggests that static correction 
modes dominate the dynamics of this example. 


9N SOLUTION 
5S SOLUTION 


— * 4N5S SOLUTION 
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COMPARISON OF COMPUTER SIMULATION TIMES 


As shown in the table below, retention of a substantial number of high- 
frequency normal vibration modes leads to very small step size and ultimately 
exceptionally large computer times. The computer times indicated are on a 
heavily loaded Prime 750 supermini computer. 


Comparison of Simulation Times 


Type 

T end 
(sec] 

CPU 

[sec] 

RMS integration 
stepsize 
[sec] 

5S solution 

2.0 

106 

0.47812E-01 

5N solution 

2.0 

401 

0.90934E-02 

9N solution 

2.0 

7471 

0.65035E-03 

4N5S solution 

2.0 

7281 

0.73491E-03 



WINDSHIELD WIPER APPLICATION 

The schematic shown is a model of an automotive windshield wiper 
assembly, in which the crank-link and two connecting links are taken as 
rigid. The left and right wiper arms are modeled as flexible bodies. The 
mechanism is driven by applying a torque to the crank link that is a function 
of motor speed. 
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FINITE-ELEMENT MODEL OF WIPER ARMS 


A modest Beam finite-element model of each wiper arm is constructed as 
shown. Friction torque, as a function of wiper tip velocity, is introduced as 
a force acting in the system, as shown. 



512 



NUMERICAL RESULTS FOR WIPER ARM TIP VELOCITY 

A flexible body solution shown in the solid line predicts vibration at a 
frequency of approximately fourteen cycles per second, relative to essentially 
the same gross motion predicted by a rigid body model of the windshield wiper 
mechanism. Experimental results with the actual system indicate an oscillation 
frequency of approximately thirteen cycles per second, very close to that 
predicted by the articulated structure model. 
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...: RIGID 
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STATUS AND DEVELOPMENTS 


The DADS flexible system dynamics code is now functioning and has been 
used to analyze a number of small and intermediate scale applications. A 
commercial version of the software is expected to be available from Computer 
Aided Design Software Incorporated, of Oakdale, Iowa. Extensions are cur- 
rently under way to enhance capability of the code to represent selected 
aspects of space structure dynamics. 


• DADS Flexible Code Is Now Functioning 

• A Commercial Version of The Code Will 
Be Available Late In 1985 

• Extensions Are Under way To Enhance 
Capabilities For Space Structure Dynamics 
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